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Abstract 

We simultaneously fit a QCD-inspired parameterization of all accelerator data on forward proton- 
proton and antiproton-proton scattering amplitudes, together with cosmic ray data (using Glauber the- 
ory), to predict proton-air and proton-proton cross sections at energies near y/s ~ 30 TeV. The p-air 
cosmic ray measurements provide a strong constraint on the inclusive particle production cross section, 
as well as greatly reducing the errors on the fit parameters — in turn, greatly reducing the errors in the 
high energy proton-proton and proton-air cross section predictions. 



*Work partially supported by Department of Energy contract DA-AC02-76-Er02289 Task B. 

tWork partially supported by Department of Energy Grant No. DE-FG02-95ER40896 and the University of Wisconsin 
Research Committee with funds granted by the Wisconsin Alumni Research Foundation. 

tWork partially supported by the U.S. Department of Energy under Grant No. DE-FG02-91ER40626. 



The energy range of cosmic ray experiments covers not only the energy of the Large Hadron CoUider 
(LHC), but extends beyond it. Cosmic ray experiments can measure the penetration in the atmosphere 
of these very high energy protons — however, extracting proton-proton cross sections from cosmic ray ob- 
servations is far from straightforward By a variety of experimental techniques, cosmic ray experiments 
map the atmospheric depth at which cosmic ray initiated showers develop. The measured quantity is the 
shower attenuation length (A,„), which is not only sensitive to the interaction length of the protons in the 
atmosphere (Ap_aii), with 
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but also depends critically on the inelasticity, which determines the rate at which the energy of the primary 
proton is dissipated into electromagnetic shower energy observed in the experiment. The latter effect is taken 
into account in Eq. by the parameter k\ nip is the proton mass and Cp^fair ^^"^ inelastic proton- air cross 
section. The departure of k from unity depends on the inclusive particle production cross section in nucleon 
and meson interactions on the light nuclear target of the atmosphere and its energy dependence. 

The extraction of the pp cross section from the cosmic ray data is a two stage process. First, one calculates 
the p-air total cross section from the inelastic cross section inferred in Eq. (Q) , where 
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Next, the Glauber method Q transforms the value of (J^pl\i^ into a proton-proton total cross section Opp] all 
the necessary steps are calculable in the theory, but depend sensitively on a knowledge of B, the slope of 



—n^ , the pp differential elastic scattering cross section, where 



B = 



dt \ dt 
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In Eq. (H) the cross section for particle production is supplemented with cr°_air ^'^'^ ^ 



the elastic and 
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quasi-elastic cross section, respectively, as calculated by the Glauber theory, to obtain the total cross section 
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We show in Fig.|| plots of B as a function of o-pp, for 5 curves of different values of o'p"fair- This 
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Figure 1: The B dependence of the pp total cross section Cpp. The five curves are lines of constant Cp'^^j^, of 440, 490, 540, 
590 and 640 mb — the central value is the published Fly's Eye value, and the others are itia- and ±2a. The dashed curve is a 
plot of our QCD-inspired fit of B against Opp. The dot is our fitted value for ^/s = 30 TeV, the Fly's Eye energy. 

summarizes the reduction procedure from the measured quantity (of Eq. 0) to Opp Q . Also plotted in 
Fig. is a curve (dashed) of B vs. app which will be discussed later. Two significant drawbacks of this 
extraction method are that one needs: 

1. a model of proton- air interactions to complete the loop between the measured attenuation length A^ 
and the cross section crp"fair' the value of k in Eq. (|l|). 
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2. a simultaneous relation between B and app at very high energies — well above the region currently 
accessed by accelerators. 

A proposal to minimize the impact of theory on these needs is the topic of this note. 

We have constructed a QCD-inspired parameterization of the forward proton-proton and proton-antiproton 
scattering amplitudes |^ which is analytic, unitary and fits all accelerator data[Q of crtot, B and p, the ratio 
of the real-to-imaginary part of the forward scattering amplitude; see Fig. 0. In addition, the high energy 




Figure 2: The simultaneous QCD-inspired fit of total cross section app, p and B vs. y/s, in GeV, for pp (squares) and pp 
(circles) accelerator data: (a) o-pp, in mb, (b) p, (c) Nuclear slope B, in GeV~^ 

cosmic ray data of Fly's Eye ||] and AGASSA experiments are also simultaneously used, i.e., k from 
Eq. (|^) is also a fitted quantity — we refer to this fit as a global fit 0. We emphasize that in the global fit, 
all 4 quantities, crtot, B, p and k, are simultaneously fitted. Because our parameterization is both unitary 
and analytic, its high energy predictions are effectively model- independent, if you require that the proton is 
asymptotically a black disk. Using vector meson dominance and the additive quark models, we find further 
support for our QCD fit — it accommodates a wealth of data on photon-proton and photon-photon interac- 
tions without the introduction of new parameters^. In particular, it also simultaneously fits app and B, 
forcing a relationship between the two. Specifically, the B vs. app prediction of our fit completes the relation 
needed (using the Glauber model) between app and Cp^^'l^jj.. The percentage error in the prediction of app at 
^/s = 30 TeV is « 1.2%, due to the statistical error in the fitting parameters (see references H, §). A major 
difference between the present result, in which we simultaneously fit the cosmic ray and accelerator data, 
and our earlier result |Q, in which only accelerator data are used, is a significant reduction (about a factor 
of 2.5) in the errors of app at = 30. 

In Fig.^, we have plotted the values of app vs. a^p_^l^^^ that are deduced from the intersections of our 
B-app curve with the cr™!^;^. curves of Fig. 1. Figure fallows the conversion of measured CTp'^i^jj. cross sections 
to app total cross sections. The percentage error in a^^^l^^^. is « 0.8 % near o-p^l^jj, = 450mb, due to the errors 
in app and B resulting from the errors in the fitting parameters. Again, the global fit gives an error of a 
factor of about 2.5 smaller than our earlier result 0, a distinct improvement. 

When we confront our predictions of the p-air cross sections (Cp'^i^jj.) as a function of energy with pub- 
lished cross section measurements of the Fly's Eye and AGASSA groups, we find that the predictions 
systematically are about one standard deviation below the published cosmic ray values. It is at this point 
important to recall Eq. (|l|) and remind ourselves that the measured experimental quantity is A^ and not 
"'p^^air- We emphasize that the extraction of a^"^]^^^. from the measurement of A^ requires knowledge of the 
parameter k. The measured depth Amax at which a shower reaches maximum development in the atmo- 
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Figure 3: A plot of the predicted total pp cross section cTpp, in mb vs. the measured p-air cross section, o"pl!?l^;j, in mb. 




Figure 4: The AGASA and Fly's Eye data for a'^^l^^, 
fit, using the common value of A; = 1.349 



in mb, as a function of the energy, y's, in GeV, as found in our global 



sphere, which is the basis of the cross section measurement in Ref . |^ , is a combined measure of the depth 
of the first interaction, which is determined by the inelastic cross section, and of the subsequent shower 
development, which has to be corrected for. Xj^ax increases logarithmically with energy with elongation rate 
(AXinax per decade of Lab energy) of 50-60 g/cm^ in calculations with QCD-inspired hadronic interaction 
models. The position of X^ax directly affects the rate of shower attenuation with atmospheric depth, which 
is the alternative procedure for extracting fJ^^fl^iy The rate of shower development and its fluctuations are 
the origin of the deviation of k from unity in Eq. (^. Its predicted values range from 1.5 for a model where 
the inclusive cross section exhibits Feynman scaling, to 1.1 for models with large scaling violations . The 
comparison between prediction and experiment is further confused by the fact that the AGASA ^ and Fly's 
Eye experiments used different values of k in the analysis of their data, i.e., AGASA used A: = 1.5 and 
Fly's Eye used fc = 1.6. 

We therefore decided to let k he a free parameter and to make a global fit to the accelerator and cosmic 
ray data, as emphasized earlier. This neglects the possibility that k may show a weak energy dependence 
over the range measured. Recently, Pryke|^ has made Monte Carlo model simulations that indicate that k is 
compatible with being energy- independent. Using an energy-independent k, we find that k = 1.349 ± 0.045, 
where the error in k is the statistical error of the global fit. By combining the results of Fig.^(a) and Fig.^, 
we can predict the variation of crp"f|^;j. with energy, y/s. In Fig.^ we have rescaled the published high energy 
data for 0'''^!^;^. (using the common value of fc — 1.349), and plotted the revised data against our prediction 
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Figure 5: A plot of the QCD-inspired fit of the total nucleon-nucleon cross section <Tyui, in mb vs. sfs^ in Gev. The cosmic 
ray data that are shown have been converted from cTp^f^ir ^pp ^^i^g results of Fig. ^and the common value of A; = 1.349, 
found from our global fit. 

The plot of Gpp vs. y/s, including the rescaled cosmic ray data is shown in Fig.|^. Clearly, we have an 
excellent fit, with good agreement between AGASA and Fly's Eye. In order to extract the cross sections' 
energy dependence from the cosmic ray data, the experimenters of course assigned energy values to their cross 
sections. Since the cosmic ray spectra vary so rapidly with energy, we must allow for systematic errors in k due 
to possible energy misassignments. At the quoted experimental energy resolutions, ALog]^Q(£'iab(ev)) = 0.12 
for AGASSA§ and ALogio(£^iab(ev)) = 0.4 for Fly's Eye§, where E^^^ is in electron volts, we find from 
the curve in Fig.| that Ak/k = 0.0084 for AGASSA§ and Ak/k ^ 0.0279 for Fly's Eye§. We estimate 
conservatively that experimental energy resolution introduces a systematic error in k such that A/cgystomatic = 
\/(AA:^QASSA + ^^flyseye)/2 = 0-028. Thus, we write our final result as fc = 1.349 ± 0.045 ± 0.028, where 
the first error is statistical and the last error is systematic. 

Recently, Pryke^ has published a comparative study of high statistics simulated air showers for pro- 
ton primaries, using four combinations of the MOCCA[0 and CORSIKAljll] program frameworks, and 
SIBYLL||l^ and QGSjet|l^ high energy hadronic interaction models. He finds k — 1.30 ± 0.04 and 
k = 1.32 ± 0.03 for the CORSIKA-QGSjet and MOCCA-Internal models, respectively, which are in ex- 
cellent agreement with our measured result, k — 1.349 ± 0.045 ± 0.028. 

Further, Pryke§ obtains k = 1.15 ± 0.03 and k = 1.16 ± 0.03 for the CORSIKA-SIBYLL and MOCCA- 
SIBYLL models, respectively, whereas the SYBILL[|l| group finds k — 1.2, which is not very different from 
the Pryke value. However, the SYBILL-based models, with k =1.15-1.20, are significantly different from 
our measurement of fc = 1.349 ± 0.045 ± 0.028. At first glance, this appears somewhat strange, since our 
model for forward scattering amplitudes and SIBYLL share the same underlying physics. The increase of the 
total cross section with energy to a black disk of soft partons is the shadow of increased particle production 
which is modeled in SYBILL by the production of (mini)-jets in QCD. The difference between the k values of 
1.15-1.20 and 1.349 results from the very rapid rise of the pp cross section in SIBYLL at the highest energies. 
This is an artifact of the fixed cutoff in transverse momentum used to compute the mini-jet production cross 
section, and is not a natural consequence of the physics in the model. There are ways to remedy this. 

In conclusion, the overall agreement between the accelerator and the cosmic ray pp cross sections with 
our QCD-inspired fit, as shown in Fig.^, is striking. We find that the accelerator and cosmic ray pp cross 
sections are readily reconcilable using a value of k = 1.349 ± 0.045 ± 0.028, which is both model independent 
and energy independent — this determination of k severely constrains any model of high energy hadronic 
interactions. We predict high energy app and 0'^'^!^;^. cross sections that are accurate to ~ 1.2% and 0.8%, 
respectively, at = 30 TeV. 

At the LHC = 14 TeV), we predict CTtot = 107.9± 1.2 mb for the total cross section, B = 19.59±0.11 
(GeV/c)~^ for the nuclear slope and p = 0.117 ± 0.001, where the quoted errors are due to the statistical 
errors of the fitting parameters. 
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In the near term, we look forward to the possibihty of repeating this analysis with the higher statistics 
of the HiRes [Q cosmic ray experiment that is currently in progress and the Auger Observatory. 
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